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Abstract: Large-mass bolometers are used in particle physics experiments to search for rare 
processes. The energy threshold of such detectors plays a critical role in their capability to search 
for dark matter interactions and rare nuclear decays. We have developed a trigger and a pulse shape 
algorithm based on the matched filter technique which, when applied to data from test bolometers 
of the CUORE experiment, lowered the energy threshold from tens of keV to the few keV region. 
The detection efficiency is in excess of 80%, and nearly all nonphysical pulses are rejected. 
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1. Introduction 

Bolometers are detectors in which the energy from particle interactions is converted into heat and 
measured via the resulting rise in temperature. They provide excellent energy resolution, though 
their response is slow compared to conventional detectors. These features make them a suit- 
able choice for experiments searching for rare processes, such as neutrinoless double beta decay 
(OvDBD) and dark matter (DM) interactions. 

The CUORE experiment will search for OvDBD of i^°Te ^ using an array of 988 Te02 
bolometers of 750g each. Operated at a temperature of about lOmK, these detectors maintain an 
energy resolution of a few keV over their energy range, extending from a few keV up to several 
MeV. The measured resolution in the region of interest (2527keV) is about SkeVFWHM; this, 
together with the low background and the high mass of the experiment, determines the sensitivity 
to the OvDBD. CUORE could also search for DM interactions, provided that the energy threshold 
is sufficiently low. DM candidates such as weak interacting massive particles (WIMPs) and axion- 
Uke particles (ALPs) are expected to produce signals at energies below ~ 30keV, and the 
interaction rate increases as the energy decreases. 

The existing CUORE energy threshold, achieved using a trigger algorithm applied to the raw 
data samples, is of the order of tens of keV. If the threshold were of few keV, the CUORE experi- 
ment could be sensitive to DM interactions and thus play an important role in this growing research 
area. At low energies the detection capability is limited by the detector noise. Electronics spikes, 
mechanical vibrations, and temperature fluctuations can produce pulses that, if not properly iden- 
tified, generate nonphysical background. The lower the energy released in the bolometer, the more 
difficult it is to discriminate between physical and nonphysical pulses. 
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In this paper we present a method to lower significantly the energy threshold, and a pulse shape 
parameter featuring high rejection power even at low energy. The advantage of this method comes 
from running the trigger and pulse shape algorithms on data that have been previously processed 
with the matched filter technique A lower energy threshold can be obtained because, compared 
to raw data, matched filtered data feature a higher signal to noise ratio. The algorithms need as 
input the expected shape of the signal and the noise power spectrum of the bolometer. No manual 
tuning is needed, because all the parameters can be set at optimal values automatically. Application 
of this technique to data from a test detector shows that the energy threshold can be lowered to a few 
keV, and that nonphysical pulses can be efficiently identified. The technique we have developed 
can be used with any kind of bolometric detector. 



2. Experimental setup 

The algorithms developed in this paper were tested on two bolometers operated by the CUORE 
collaboration at the LNGS underground laboratory. The main purpose of the test was to check one 
of the first production batches of CUORE crystals [Q]. In the following we will describe our studies 
taking the first bolometer as example, while the final results will be shown for both. 

A CUORE bolometer is composed of two main parts, a Te02 crystal and a neutron transmu- 
tation doped Germanium (NTD-Ge) thermistor [^, ^]. The crystal is cube-shaped (5x5x5 cm^) and 
held by Teflon supports in copper frames. The frames are connected to the mixing chamber of 
a dilution refrigerator, which keeps the system at the temperature of 10 mK. The thermistor is 
glued to the crystal and acts as thermometer. When energy is released in the crystal, its temperature 
increases and changes the thermistor's resistance [^, ^. The thermistor is biased with a constant 
current, and the voltage across it constitutes the signal [[To|]. Usually, a Joule heater is also glued 
to the crystal. It is used to inject controlled amounts of energy into the crystal, to emulate signals 



produced by particles [|ri|, |12[ ]. 

The typical response of CUORE bolometers to particles impinging on the crystal is of order 
100;UV/MeV. The signal frequency bandwidth is — 20 Hz, while the noise components extend 
to higher frequencies. The signal is amplified, filtered with a 6-pole active Bessel filter with a 
cut-off frequency of 12 Hz and 20 Hz for the first and the second bolometer, respectively, and then 
acquired with an 18-bit ADC with a sampling frequency of 125 Hz. Given the low event rate and 
the small frequency bandwidth of the signal, it is possible to save the complete continuous data 
stream of the detector voltages, so that it can be reprocessed offline. An online software trigger 
is used as well: when it fires, a fixed number of data samples is saved. There are three kinds of 
triggers: those caused by particle-like pulses, random waveform samplings, and those generated 
by the Joule heater. Each event is tagged with a flag that indicates the trigger type. The standard 
online particle trigger algorithm is based on a simple derivative approach: it fires when the slope of 
the signal remains above threshold for a certain amount of time. In the measurements discussed in 
this paper, the energy threshold reached with this algorithm was 1 1 keV on the first bolometer and 
40keV on the second one. 

A typical signal produced by an 88keV 7 particle generated by a '^Vmrpg de-excitation in 
the crystal is shown in Fig. |l|. The length of the acquisition window was 5.008 s, correspond- 
ing to M = 626 samples. The energy calibration was made by exposing the bolometer to a ^^^Th 
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Figure 1. Signal produced by a 88keV / particle fully absorbed in the detector. The signal is sampled at 
125 Hz and the length of the window corresponds to 5.008 s (626 samples). 
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Figure 2. Energy spectrum of the calibration measurement. All lines are generated by the ^^-^Th calibration 
source except for the lines at 1461 keV, due to '*'^K contamination of the refrigerator, and at 5407 keV, due 



to ^'"Po contamination in the Te02 crystal. 



source; for what concerns our studies, the calibration function can be considered linear at ener- 
gies below lOOkeV. For the first bolometer the conversion factor from ADC counts to energy 
is 15 ADCcounts/keV. The acquired calibration spectrum is shown in Fig. ^ In it, the most 
prominent peaks of the ^■'^Th source are visible, together with two lines at 1461 keV and 5407 keV, 
which are due to '^''K contamination of the refrigerator and ^'"Po contamination of the crystal bulk, 
respectively. The ^^'^Po line decays with a half-life of 138 days and is visible in recently-produced 
crystals. In the data analyzed in this paper, the rate of the ^'"Po line was about lOmHz and signifi- 
cantly affected the detection efficiencies at low energies, as will be shown in the next sections. 
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3. Data filtering 



The trigger and the pulse shape identification algorithms presented in this paper operate on data 
samples processed using the matched filter technique This filter is designed to estimate the 
amplitude of a signal, maximizing the signal to noise ratio. The transfer function is matched to 
the signal shape, and pulses with different shape are suppressed. The noise power spectrum of the 
detector, N{(Ok), and the signal shape, Si, are needed to build the transfer function: 



(3.1) 



where s{(Ok) is the Discrete Fourier Transform (DFT) of Si, im is the maximum position of 5, in the 
acquisition window, and his a normalization constant that leaves unmodified the amplitude of the 
signal: 



The noise power spectrum at the filter output can be estimated as 

Nf{cOk)=h^ 



from which the resolution Of can be derived: 



o} = Y^Nf{(Ok)=h. 

k 



(3.2) 



(3.3) 



(3.4) 



The estimation of N{(Ok) is made by averaging the power spectrum of a large set of data windows 
not containing signals, while the signal shape si is obtained by averaging a set of particle signals 
selected from the calibration measurement in the 1-3 MeV range. The filtered average pulse is 
symmetric and has the same amplitude as the original (see Fig. ||, left). On the other hand, noise 
frequencies with low signal to noise ratio are suppressed (see Fig. ^ right). The resolution before 



the filter, a, is 1.1 keV and, according to Eq. [3.4[ , is reduced by the filter to a/ = 0.3 keV. 
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Figure 3. Average pulse (left) and noise power spectrum (right) of the detector, before (solid black line) and 
after (dashed red line) the filter. The filtered average pulse is symmetric and has the same amplitude as the 
original. The filtered noise is enhanced in the signal frequency bandwidth and suppressed elsewhere. 
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The filtering of the data samples can be implemented in two ways, in the time domain or in 
the frequency domain. In the time domain one has to implement a finite impulse response (FIR) 



filter whose kernel is the DFT of Eq. |3JJ. In the frequency domain, data samples first have to be 
transformed using DFT, then multiplied by Eq. [3.1[ , and then transformed back in the time domain. 
The number of floating-point operations of the two methods is of the same order of magnitude, 
namely 0{M^). However, in the frequency domain convolution, EFT algorithms can be used in 
place of the standard DFT, reducing the number of operations from 0{M^) to O(MlogM). In view 
of a large scale application, we chose to filter the data in the frequency domain. 

DFT algorithms assume that the window to be filtered is periodic. If the input data are not 
periodic, the convolution will falsely pollute the left side of the window with data from the right 



side and vice versa. This problem is known as "DFT convolution wraparound" [13] and can be 
solved by padding the response function of the filter with zeros, such that a slice of the window is 
correctly filtered. We follow these steps to build the transfer function: 

1. Compute the matched filter transfer function H{(Ok) of length M using the average signal 
shape Sj and the noise power spectrum N{(Ojf:), both of length M. 

2. Transform H{(Ok) into the time domain, obtaining the response function Hi of length M. 

3. Insert M zeros in the middle of Hi, obtaining the response function Hf of length 2M. 

4. Smooth Hf in proximity of the zero insertion, obtaining the smoothed response function Hf\ 

5. Transform Hf^ into the frequency domain, obtaining the final transfer function H'^'^{(Ok), of 
length 2M. 



The smoothing of Hf is needed to avoid the Gibbs phenomenon [|14p : if Hf does not approach zero 



id . 



where the zeros will be inserted, a discontinuity will be created, introducing fake oscillations in 
its DFT. The smoothing is made by multiplying Hf by a quarter-period cosine shape. The cosine 
period is L = 120 samples, and the quarter-period that goes from one to zero is used to smooth 
Hf at the left of the zero insertion, while the quarter-period that goes from zero to one is used to 
smooth Hf at the right of the zero insertion. 

When filtering a data sample of length 2M, the first and the last M/2 samples will be spoiled 
while the M middle samples will be correctly filtered (Fig. 0). The continuous data flow is pro- 
cessed by filtering windows of length 2M overlapped by M and taking the M middle samples of 
each window as output. The concatenation of the middle samples produces a continuous flow of 
correctly filtered data (Fig. ^). 

4. Trigger algorithm 

Some slices of the filtered data are shown in Fig. |6| and exhibit the following features: 

1 . The noise fluctuations are reduced. 

2. The baseline is zero. 
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Figure 4. Padding of the filter response function. To avoid the DFT wraparound problem the original 
response function is padded in the middle with a number of zeros equal to its length. In this way the product 
of the response function DFT and of the data samples DFT produces a correct convolution in correspondence 
of the response function zeros. 
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Figure 5. Filtering of the continuous data flow. Each window has size 2M (white rectangle), and the M 
middle samples are correctly filtered (gray rectangle). The windows are overlapped by M, such that the M 
middle samples of each window are concatenated. 



3. The filter is sensitive to the shape of the expected signal and pulses with different shape are 
suppressed. 

Pulses are triggered when the data samples exceed a positive threshold, and a new trigger is 
possible when the data samples return below threshold, or after that a local maximum is found. 
The trigger threshold d is defined in terms of number of sigma of the filtered noise {Of), which is 
known a priori (see Eq. 3^). The triangles in Fig. ^ show the found triggers and their positions, 
using the threshold d = 4.5a/. The trigger position is then shifted with respect to the maximum, so 
that it lies in the middle of the rise of the pulse. This operation is needed since in the data analysis 
the trigger is expected to lie before the maximum of a signal. The size of the shift back is fixed 
and is computed from the time difference between the middle of the rise and the maximum of the 
non-filtered average pulse. 

Although the matched filter dramatically simplifies the implementation of trigger algorithms, 
there is a drawback. When high-energy pulses are triggered, a set of secondary pulses is seen. This 
is due to the fact that the filtered pulse has two symmetric lobes above zero, that would be marked 
as signal by the trigger (see Fig. This problem is solved as follows. The amplitude above zero 
of the side lobes is a constant fraction a of the amplitude of the main lobe, and is estimated from 
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Figure 6. Windows of data samples containing pulses. The raw data (solid black line) have been shifted such 
that the first sample has zero ADC counts; the filtered data (dashed red line) are not modified. Red triangles 
identify the pulses detected by the trigger algorithm. The filter removes the baseline drifts and suppresses 
pulses with different shapes than the expected signal. 
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Figure 7. Filtering of a high-energy signal (5407 keV a). When a high energy pulse occurs (left) the side 
lobes of the filtered signal exceed the threshold, and fake signals could be triggered. To avoid this, the 
regions corresponding to the lobes are vetoed and no fake trigger is saved. In the right picture: the left lobe 
of the pulse, and the vetoed region (gray band). 



the filtered average pulse. When a signal with amplitude A > 0/a is triggered, the data samples 
corresponding to the expected positions of the side lobes are vetoed, and no trigger is allowed 
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in those regions (see Fig. This is an important source of inefficiency for the measurements 
discussed in this paper, since the ^^°Po contaminations were high and the amount of dead time was 
dominated by the rate of the 5407 keV line. When crystals are sufficiently aged (one or two years 
since production), the rate is so low that the dead time is neghgible. 

5. Detection efficiency 

To measure the detection efficiency, an energy scan with the external heater was performed using 
sequences of pulses from the keV region up to ~ 100 keV (Fig. ^). 
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Figure 8. Energy spectrum of heater pulses in the scan measurement. Only pulses flagged as "heater" by 
the DAQ are shown. Below ^ 3keV, pulses were fired at 3 different energies. 

Comparing the trigger output with the DAQ heater flag, each pulse is considered detected if: 

1 . The trigger fired on a data sample in a time window that goes from 20 samples before the 
flag to 30 samples after. 

2. The energy obtained from the matched filter amplitude is compatible with the energy ex- 
pected for the specific heater pulse. 

In order to verify this last condition, energy spectra of events related to each heater amplitude set 
are considered. In each spectrum there is a single peak and some outlier in which the matched 
filter failed to assign the correct amplitude, due to noise fluctuations or pile-up phenomena. Using 



the resolution Of in Eq. 3.4, all the events not included in a 2FWHM energy window (4.7a/), 
centered on the mean energy of the peak, are flagged as not detected. The detection efficiency is 
then estimated as the ratio of the events within this energy window and the total number of heater 
pulses fired at that energy. 

The measured efficiency, compared with the one of the standard trigger, is shown in Fig. ^ 
for the two bolometers. The efficiency increases with energy and reaches a plateau. In the first 
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Figure 9. Detection efficiency computed via the heater scan for the standard trigger (open black circles) 
and the new trigger (solid red circles). The efficiency increases with energy, reaching a plateau. In the first 
bolometer the plateau is reached at 1 1 keV by the standard trigger and at ^ 2.5 keV by the new trigger. In 
the second bolometer the plateau is reached at 30-^50keV by the standard trigger and at ^ 5keV by the 
new trigger. 



bolometer the plateau efficiency amounts to 90% and is reaclied at ~ 2.5keV. Tiiis 10% efficiency 
loss is compatible with the veto dead time due to the ^'"Po line rate. The same efficiency plateau 
is reached by the standard trigger at ~ 1 1 keV. In the second bolometer the improvement is even 
more evident: the 84% plateau is reached at 30-^50keV by the standard trigger and at ~ 5keV 
by the new trigger. 



6. Pulse shape identification 

As previously stated, the matched filter is sensitive to pulses having the same shape as the expected 
signal. When a pulse with different shape occurs, its filtered amplitude is suppressed and the filtered 
shape also is different from the expectation (see Fig. To suppress fake signals we implemented 
a pulse shape indicator. This indicator is not used to determine the presence of a pulse, but can be 
used in the data analysis to reduce the nonphysical background. 

The algorithm is relatively simple. A triggered pulse is fitted using a cubic spline of the filtered 
average pulse, and the ;^^/ndf of the fit is used as the shape indicator. The cubic spline is necessary 
to produce a continuous function and to fit fractional time delays. The fitted parameters are the 
amplitude and the position of the pulse, while the baseline is fixed at zero. Since the matched 
filter is itself a fit of the data samples, this fit only serves to remove the digitization effects and 
to estimate the x^- Before fitting, the maxima of the fit function and of the triggered pulse are 
aligned. The pulse position is then varied in the range [—1,1] samples, which is the maximum shift 
due to the digitization. The amplitude is varied in the range [A,A(1 +£)], where A is the amplitude 
of the filtered signal and e is the maximum allowed spread in amplitude due to the digitization. 
This parameter is estimated from the filtered average pulse, computing the difference between the 
amplitude of the maximum (that is equal to one) and the amplitude of the consecutive sample. Some 
fits to pulses triggered on the first bolometer are shown in Fig. [l^ The fit range L corresponds to 
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Figure 10. Fits of pulses triggered on the first bolometer; points represent data and residuals, solid lines 
represent the fit functions. Left: an event in the 3-4 keV region, SI = 0.5. Right: another event in the 3-4 keV 
region, SI = 23. 



four FWHM of the filtered average pulse. Finally, the shape indicator (SI) is computed as 

SI = Vi^Z:^ (6.1) 

where y-^ is the filtered signal, / is the minimized fit function, and Ol is the amount of noise 
expected in a window of length L. The error of each point, in fact, is not the integral error of the 
entire filtered window in Eq. 3.4, but it is smaller since low frequencies are not seen in a short 
window: 

It must be remarked that, even though Eq. il has the form of a it does not follow a true 
distribution. Although the expected value is still 1, the variance is not 2/(L — 2) since the noise 
fluctuations are correlated at the filter output (see residuals in Fig. 10, left). 

The distribution of SI versus energy from 4.6 days of data acquired without calibration sources 
is shown for both bolometers in Figs. [Tll-|l2[ The signal region is well separated, and nearly all 
the non-physical pulses can be removed with a cut on this variable. We do not discuss here the 
optimization of this cut, which depends on a specific analysis. 



7. Conclusions 

We present a method to lower significantly the energy threshold of bolometric detectors, and a 
pulse shape parameter featuring high rejection power even at low energy. This result is achieved 
by running the trigger and pulse shape algorithms on data that have been previously processed with 
the matched filter technique. All the input parameters required by these algorithms can be set at 
optimal values automatically and no manual tuning is needed. 



-10- 



10 



: Bolometer 1 


... •*. 












signal 

' .1. 




. , • ■ ■ 



10 



Energy [keV] 



10^ 



lO'*! 



C/5 



DUiUlllCLCl Z. 














signal 











10 



Energy [keV] 



10^ 



Figure 11. Shape indicator as a function of energy for the first and the second bolometer The band at low 
SI is populated by signal events, where the lines at 30keV (X-ray from Te) and at 88keV (7 from '^^'"Te) 
are visible. At higher values there are triggered mechanical vibrations and spikes. 
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Figure 12. Shape indicator distribution for the first and the second bolometer. The distribution of all the 
events in the range 2.5-lOOkeV for the first bolometer, and in the range 5-100 keV for the second bolometer 
(black solid line), shows a good separation of the signal region (peak around logio(Sl) = 0), which is even 
more evident if only the events in the range 10-lOOkeV are selected (red dashed line). 



The application to two Te02 bolometers lowered their thresholds from 1 1 and 40keV to 2.5 
and 5keV, respectively, with a detection efficiency in excess of 80%. The pulse shape algorithm we 
developed allows the rejection of nonphysical pulses, such as spikes and temperature fluctuations. 
With these tools, the TeOi bolometers of the CUORE experiment will be able to search for new 
physics at low energies, such as rare nuclear decays and dark matter interactions. 
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